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to use formulae from continuum ChPT. Finally, by dropping the data points near 
135 MeV from the fits, we can explore the range of pion masses that is needed in 
SU(2) ChPT to reliably extrapolate to the physical point. 
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I. INTRODUCTION 

One of the most fascinating aspects of QCD [lj], the theory of strong interactions, is the 
non- analytic behavior of its Green's functions close to the chiral limit, that is with two or 
three quark masses taken small, m q <C A 2 /B, where A~ 1 GeV is a typical hadronic scale 
and B is a condensate parameter which we will define (and determine) below. 

The structure of these non-analytic contributions can be worked out in the effective field 
;heory approach which, for the case at hand, is known as chiral perturbation theory (ChPT) 
. In this setup physical quantities are expanded in powers of p/A, where p is an external 
momentum, and the quark mass is treated as m q ~ 0(p 2 ). Depending on whether this is 
done for the two {u, d) or three {u, d, s) lightest flavors, the framework is known as SU(2) or 
SU(3) ChPT. 

In either case, at the leading order of the chiral expansion there are two low-energy 
constants, defined from the pion decay constant as 

F= lim F w , B= lim {-{0\qq\0)/F w } (1) 

where q denotes one of the light flavors (i.e. q = u or q = d) in the SU(2) case, and similarly 
(with also m s sent to 0) in the SU(3) case. At the next-to- leading (NLO) order, i.e. at 0{p A ) 
in the chiral counting, 7 new low-energy constants show up in the SU(2) framework, or 10 
low-energy plus 2 high-energy constants in the SU(3) framework. These low-energy constants 
parametrize the above mentioned chiral logarithms in the Green's functions of QCD. Their 
numerical values can be determined either from experiment or from an ab-initio solution of 
QCD in the relevant (small coupling and light quark mass) regime, as is provided by lattice 



QCD J4-6]. Since the chiral logarithms show up as rather subtle effects, meaningful results 
can only be obtained from lattice data which have excellent statistical precision and explore, 
at the same time, a wide enough range of lattice spacings, quark masses and box volumes 
such that all sources of systematic error can be controlled and eventually removed. 

In this paper we provide such a determination of the SU(2) low-energy constants £3 
and £4. Their numerical values are extracted from the quark mass dependence of M T and 
F n , respectively, and complemented by numerical values of the leading order low-energy 
constants F and B. We use staggered fermion simulations with Nf = 2 + 1 dynamical 
flavors, that is two degenerate light quarks of variable mass mi, and an active strange quark 
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whose mass m s is pinned down at its physical value. As a result of this our values of the 
low-energy constants are supposed to coincide with those in the real word. A preliminary 
account of our work (based on a smaller dataset) was given in j^]. 

The remainder of this paper is structured as follows. In Section [XT] we specify the gauge 
and fermion actions used and list the ensembles which go into the determination of the chiral 
low-energy constants. Furthermore, details are given how we calculate the pion mass and 
decay constant, and how we correct the latter for the effect of the finite spatial volume of the 
box (which is always a small correction, since our data satisfy ?>.'5<M. K L<Q.&). In Section 
DTD we specify the procedure through which we set, for each bare coupling (3 = 6/<7o, the 
lattice spacing a and the physical values of the bare quark masses m; = (m M +md)/2 and m s . 
Section IIVI contains the core part of the present investigation, an analysis of our data with 
SU(2) ChPT at NLO, with details of how we select adequate mass windows and determine 
the systematic uncertainty of the fitted low-energy constants. Section |V] contains a similar 
though less mature analysis at NNLO, where again the main goal is to determine the NLO 
coefficients, with and without the help of some priors on the remaining NLO and NNLO 
low-energy constants. This helps to give a reliable estimate of the theoretical uncertainty of 
the NLO results obtained. A summary and a comparison with the findings of other recent 
lattice studies of SU(2) NLO low-energy constants is presented in Section [VTl 



II. LATTICE DATA 

In this section we specify the lattice actions used, list the ensembles which go into the 
determination of the SU(2) low-energy constants, give details of how we extract the pion 
mass and decay constant on a given ensemble, and describe the procedure by which we 
remove the (small) impact of the finite spatial box size L on the data. 



A. Lattice action and ensemble generation 



The lattices are generated with a tree-level Symanzik improved gauge action [8| and 2+1 
flavors of staggered quarks with 2 levels of stout-smearing [9|. The action is specified in full 



detail in Ref. 



10|. The a 



improvements (see e.g. 



gorithm used is a combination of HMC and RHMC with standard 



for an overview). Some ensembles were generated for scale 
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setting purposes in previous finite temperature studies [10Ul2Hl6j]. and some were generated 
specifically for the present investigation. 

We adopt a mass independent scale setting, that is the lattice spacing a depends only on 
the coupling in the gauge action, not on the quark masses mi,m s . With this choice it is 
straightforward to adjust, for each (3, the strange quark mass roughly to its physical value by 
tuning the ratio Mgs/M^ to its physical value, where Mg s = 2Mf C — M^ and is measured 
from the connected contribution. Starting from the symmetric point mi = m s = mf ys , one 
can lower the light quark mass mi, at fixed m s , until the ratio M^/Mgs assumes its physical 
value. This is one possible definition of the physical point; other definitions differ from this 
one just in the choice of which ratios are affected by cut-off effects and which are not (the 
latter are those which are used in the definition of the physical point). 

From a more practical point of view it suffices to say that we simulate, for each /3, a 
number of (m;, m s ) combinations, where m s is held fixed and is close to whichever definition 
of the physical strange quark mass that one may adopt, while m/ varies between roughly the 
physical light quark mass and 4 to 10 times this value (depending on f3). The precise value 
of mf hyh is determined, a posteriori, by means of an interpolation, as described in Section 
IIHI below. A summary of our ensembles and their bare parameters is given in Table HI 



B. Calculating meson masses and decay constants 

A specific advantage of staggered fermions (or of any other discretization with some form 
of chiral symmetry) is that the decay constant / = \/2F of a pseudoscalar meson (in the 
following we will distinguish the two normalizations by using either the upper-case or the 
lower-case symbol) to the zero component of an axial current can be extracted without 
recurrence to any lattice-to-continuum matching factor. 

We start from the two-point function Cpp(t) between two point-like pseudoscalar density 
operators (at least one of which is projected to zero spatial momentum) which, for an 
intermediate window of the Euclidean time t, takes the form 

C PP {t) = A PP f exp(-Mf) + exp(-M(T-t))l (2) 



with T the lattice extent in the fourth direction. The mass M corresponds to the mass of 
the lightest asymptotic state with the right quantum numbers (here ir or K), whereas the 



5 



p L 3 x T 



ami 



am^ 



(3 L 3 x T ami 



am q 



3.45 24 3 x 32 
16 3 x 32 
12 3 x 28 
12 3 x 28 
12 3 x 28 



0.0057619 0.1573 

0.0172857 

0.0288095 

0.0403333 

0.0518571 



3.55 24 3 x 32 
16 3 x 32 
16 3 x 32 
12 3 x 28 
12 3 x 28 



0.00374878 0.1023417 

0.01312073 

0.01874390 

0.02624146 

0.03373902 



3.67 32 3 x 48 
24 3 x 32 
16 3 x 32 
16 3 x 32 
14 3 x 32 



0.00231904 0.06330976 

0.00927616 

0.01391424 

0.01739280 

0.02203088 



3.75 48 3 
40 3 
32 3 
40 3 
32 3 



x 96 0.00172000 0.048 
x 64 0.00240000 
x 64 0.00342857 
x 64 0.00480000 
x 64 0.00685000 



3.792 48 3 
40 3 
40 3 
40 3 
32 3 



x 64 0.00160714 0.045 
x 64 0.00225000 
x 64 0.00321429 
x 64 0.00450000 
x 64 0.00674300 



3.85 48 3 x 64 0.00144606 0.0394774 
40 3 x 64 0.00197387 
40 3 x 64 0.00281981 
48 3 x 64 0.00394774 
32 3 x 48 0.00578424 
32 3 x 48 0.00867636 
24 3 x 48 0.01156848 
24 3 x 48 0.01446060 



TABLE I. Overview of the staggered 2+1 flavor ensembles used in this work. 



amplitude A PP is proportional to the squared matrix element, i.e. A PP oc (0\P\xy) 2 /M, 
where \xy) denotes the pseudoscalar state put together from flavors x, y (here x,y — l,s). 

In practice it means that we determine, in a first step, the mass and the amplitude from 
the PP correlator. We do this either via the effective mass and amplitude method where 
the determination of 

1 (C(t-1) + \/C(t-l) 2 - C(T/2) 2 \ 
M eS (t) = log -i ' V ( 3 ) 

Acs{t) = exp[-M cff (t) t] + eip[-M cff (t) (T-t)] (4) 
is followed by a fit to a constant over some time region t G [t m \-n '■ t max \, or using a direct fit 
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of the correlator to the functional form (|2"1) . In either case the data are symmetrized about 
T/2, and t max <T/2. 

The decay constant is defined as f xy = (0\A4\xy) /M xy , and via the PCAC relation this 
is transformed into 

A0\P\xy) 

fxy = (m x +m y ) (5) 

xy 

where m liS denotes the quark mass of the flavor x or y. Putting things together, it follows 
that the decay constant may be obtained from the amplitude and the mass as 



f xy oc (m x +m y 



M 3 y ' 

xy 



where the missing prefactors (e.g. L 3 ) reflect normalization conventions which depend on 
the geometry, but not on the quark masses. 

For the interim step, i.e. the determination of M e fj, A e d from the correlators, a typical 
plateau is shown in Fig.[TJ A careful look reveals some zig-zag of the data close to the 



mid-time point T/2. This is commonly attributed to a back-propagating parity partner 17] 
and reflects an effect which is specific to the staggered discretization. An advantage of the 
symmetric definition of the effective mass ([3]) is that M e g(t) for an odd value of t uses only 
data from the original correlator at even t and vice versa. Accordingly, we can compare 
the results of (i) a plateau average of M eS (t) for odd t, (ii) a plateau average of M e g{t) 
for even t, and (Hi) the result of a direct fit to the ansatz (T5]) [which does not distinguish 
between even and odd timeslices]. We have carefully analyzed the impact of these options 
and found them completely insignificant compared to both the statistical uncertainty and 
(even more so) the theoretical uncertainty inherent in the precise choice of the masses and 
lattice spacings included in the chiral fit. The latter represent relevant options that will be 
discussed in detail in Sections [TV] and [V] below. 



C. Finite volume corrections 

Pseudoscalar masses and decay constants experience a systematic shift due to the finite 
spatial box length L. Approximate 3-loop and 2-loop expressions have been given for the 



ratios M 7r (L)/M^(oo) and / 7r (L)// 7r (oo), respectively, in Ref. [18|. For the range of quark 



masses used in this work they are supposed to give a reliable estimate of these (small) shifts. 
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FIG. 1. Effective masses (left panel) and amplitudes (right panel) from the pion correlator at 
f3 = 3.75, ami = 0.00172, am s = 0.048, including the plateau-fits (blue lines). 



The numerical values thus obtained vary between 0.1 and 2.7 per-mille for the pion mass, 
and between 0.2 and 7.5 per-mille for the pion decay constant. For each ensemble we thus 
calculate these ratios and apply them to our data. In the following, only the finite-volume 
corrected data are used. 



III. DETERMINING THE LATTICE SCALE AND PHYSICAL QUARK MASS 

As indicated in the previous section, we wish to determine for each j3 the lattice spacing 
a and the physical value of ami. Since the simulation points for one /3 are at a fixed value 
of am s (which is tuned to its physical value, see Sec. lII A|) . it is clear that the observables 
to be used shall include only the light but not the strange flavor, and the obvious choice is 
thus = y/2F n and M v . 

On a more technical level we proceed by means of a two step procedure. First, we ex- 
trapolate the ratio (aMu) 2 / (afu) 2 = (aM 7T ) 2 /(af n ) 2 of the squared pion mass and decay 
constant to its physical value, 1.06846. Here we use the isospin averaged and electromag- 
netically corrected pion mass M% hys = 134.8 MeV from FLAG [l9], and the PDG value of 



the pion decay constant /£ hys = 130.41 MeV 20]. In this step the purpose of the square is 



to reduce the amount of curvature, and we interpolate the data by means of a low-order 
polynomial and rational ansatz (typically with 3 parameters applied to the 5 lightest data 
points, i.e. with 2 degrees of freedom). We stress that the point where this ratio assumes 
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the desired value, 1.06846, is always very close to the lightest simulated quark mass. In 
view of this it should not come as a surprise that the values of amf hys that stem from the 
polynomial and the rational fit are always very close to each other (on the scale set by the 
statistical error). We use the average of the two as our central value; the difference should 
be seen as indicative of the systematic uncertainty of amf hys from this set of ensembles. 

In the second step we consider af^ as a function of ami. Again, we interpolate the data 
with the same polynomial and rational ansatz, and determine the ordinate value at the 
abscissa point that was specified in the previous step. This value af n is then identified with 
the product of the lattice spacing a and the PDG value f% hys = 130.41 MeV j20|; this yields 
the lattice spacing a in fm for the lattice theory at that particular value of the coupling 
0. A typical example of this two-step procedure is shown in Fig.[2j in the relevant range 
(close to the lightest mass point) the difference between the two ansatze is invisibly small. 
Furthermore, in Tab. [TT] the results for the physical light quark masses and lattice spacings 
obtained by this method are displayed. 

As a final comment, let us remark that already in these two steps one could, in principle, 
use ChPT. We rather prefer to stay with the simple yet robust procedure as sketched above. 
It makes it apparent that some of the chiral fits (with an inappropriate fitting window) going 
wild (as will be discussed in Sections IIVI and |V] below) is not linked to a potential mishap 
in the physical mass and scale determination. In other words, we take the lattice spacing 
and the physical light quark mass from an "ideal" simulation where M 7T /f n is exactly tuned 
to its physical value, and use this knowledge regardless of how many data-points enter the 
chiral fits described below. Needless to comment that for those NLO and NNLO fits which 
work fine (and which include the lightest data points), their intrinsic physical mass and 
scale determination was always found to be in very good agreement with the result of the 
procedure described in the previous two paragraphs. 



IV. FITS TO NLO-SU(2) CHPT 

In this section we will describe how we fit the mass dependence of the meson decay 
constant and its squared mass to the prediction of next-to-leading order (NLO) ChPT and 
in this way obtain the low-energy constants (LECs) appearing in these ChPT-formulae. 
Since in the simulations considered here the strange quark mass was fixed to its physical 
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FIG. 2. Fits to the ratio (aM 7r ) 2 /(a/ 7r ) 2 (left panel) and the decay constant af % (right panel) 
at P = 3.85 to obtain the physical light quark mass am phys and the lattice scale 1/a. Shown are 
quadratic (solid blue lines) and rational (dashed red lines) fits including the black data points, green 
data points were excluded from the fits. Error bands are shown by dotted lines. The dashed grey 
lines mark the value of the ratio at the physical point and the physical light quark mass extracted 
therefrom. 



p 


am phy 


s 


1/a/GeV 


a / fm 


3.45 


5.771(18) • 


10~ 3 


0.69468(67) 0.28406(27) 


3.55 


3.612(18) • 


10" 3 


0.9165(12) 


0.21531(28) 


3.67 


2.191(15) • 


10" 3 


1.3063(25) 


0.15105(29) 


3.75 


1.6889(55) 


• 10~ 3 


1.6288(15) 


0.12115(11) 


3.792 


1.5355(70) 


• 10~ 3 


1.7935(22) 


0.11002(14) 


3.85 


1.3430(98) 


• 10~ 3 


2.0410(35) 


0.09668(16) 



TABLE II. Physical light quark mass om phys , lattice scale 1/a and spacing a at different ft obtained 
by the method described in Sec. IIII1 
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value, we will only deal with the light quark mass dependence described by SU(2) ChPT. 
The LECs extracted in this way therefore contain the correct contribution of the effects due 
to the strange quark present in nature. Further, we are restricting ourselves to the case of a 
degenerate light quark mass, which we will denote by m (or am in lattice units) throughout 
the remainder of this paper, i.e., from now on we drop the subscript I on the light quark 
mass parameter. At the physical point, this mass corresponds to the average mass of the 
two light quarks observed in nature: m phys = (mj h),s + mJ^ hys )/2. 

In the following we will try to fit our data to continuum ChPT and not consider variants 
of ChPT, which take into account lattice discretization effects and/or taste violations present 

finn 

in the staggered formulation, see e.g. [171 121L |22j. Whether a continuum ansatz is suitable 



to describe our data is not a priori clear and needs to be tested. It is valid if, within the 
statistical precision of our data, no cut-off dependence is seen in the observables considered. 
As will be shown below, this is indeed the case for our M n and f n data in the relevant region. 
This can be traced back to the combined effect of (a) the two levels of stout smearing in 
the action (see Sec. HI]) and (b) the specific choice of our scaling trajectory, i.e. the scale and 
the physical quark mass being set through f n and M n (see Sec. lIIII) . While (b) ensures that 
discretization effects on both M n and f n vanish at the physical point, (a) keeps them small 
in its vicinity due to the suppression of taste violations (see e.g. 23, 24] and [ill). As we 
will see below, the combination of these two effects leads to a suppression of discretization 
effects within our statistical precision over a sufficiently wide region of quark masses and 
lattice spacings. 



A. Methods and fit formulae 



ChPT up to NLO predicts for the decay constant f n and the squared mass M% of a meson 
consisting of two mass-degenerate quarks of mass m the following functional form 

X , X 



Mi 



X 



1 + 



/* = / 



16vr 2 / 2 A3 
8^P g A| 



2. 3|, [U 
(7) 



X 



2Bm. 



(8) 
(9) 



At this order four LECs appear: the decay constant / in the SU(2)-chiral limit (m — > 0), the 
mass-parameter B and two low-energy scales A 3 and A 4 . The mass parameter B depends like 
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the quark mass itself on the renormalization scheme but the combination x of this parameter 
and the quark mass are renormalization scheme- independent and it is this combination which 
will be used exclusively in this work. The low-energy scales are related to the LECs £3 and 
£4 at the scale of the physical pion mass M phys via 

A? 



log 



3,4 



(10) 



(M? hys ) 2 ' 

Since we used the physical values / phys and (M phys ) 2 to set the scale and determine 
the physical light quark mass m phys for each set of lattice ensembles with a given gauge 
coupling (3, our ansatz should reproduce the physical point. Note also, that each set of lattice 
ensembles contained one simulated point in close vicinity of the physical point. Therefore, 
we could impose the following constraints to the chiral formulae (Eqs. (J7J |8])) 



Mi 



iPhya 



( M Phy: 



s\2 



iphys 



fP 
J 7T 



phys 



(11) 



As it is easily derived, these two constraints result in the relations 



log 
log 



x phys 


167T 2 / 2 


AI 


( x phys)2 


x phys 




AI 


xphys V 



((M phys ) 2 

f _ f phys 



X 



phys 



(12) 
(13) 



between the LECs, where ^ phys = 2Bm phys . Using them to eliminate, e.g., the low-energy 
scales from the NLO ChPT formulae, the constrained formulae can be written as 



Mi 



X 



1 + 



1 - 



X 



16n 2 f 2 

X 



log- 



X 



phys 



8vr 2 / 2 



log- 



X 
X 



+ 



X (Mf ys ) 2 - X 



phys 



X 

X 



phys 



X 



phys 



X 



phys 



X 



phys 



f _ jphys 



(14) 

(15) 



Note that x/x phys = fn/m phys . These formulae now depend on two LECs: B and /, and the 
two physical input values M% hys and / phys . In our fitting procedure only ^ phys and / will be 
treated as free parameters. For that reason, we also like to refer to these fits as parameter 
reduced fits. (Of course, by treating M phys and / phys as free parameters as well, one would 
recover the unconstrained Eqs. (JTJ, [8j) , respectively.) 

In the following, we want to perform combined (i.e., fitting M 2 and f n simultaneously) 
global fits to our lattice data at the different gauge couplings available to us. For this reason 
we make use of the lattice scale and physical light quark mass determined beforehand (see 
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Sec. lIIIj) . The meson mass and decay constant in lattice units at a given gauge coupling (3 
are converted into physical units by 

Ml = (l/^Kj, U = {l/a) p (afjp, 

respectively. Furthermore, we rewrite the combination \ = 2Bm as 

and will determine only the renormalization-independent combination ^ ph y s = 2Bm phys in 
our fits. The factor m phys will be removed, based on external data, in Sec. lVIl 



B. Combined global fits 

I. NLO-ChPT without constraints 



We begin the discussion of the chiral fits with the results of applying the unconstrained 
fit formulae, Eqs. ([71 EJ), to our data. The combined fit has four free parameters: ^; phys = 
2Bm phys , f, A3, and A4. A priori, it is not clear whether all simulated lattice spacings will 
lie in the scaling region, especially for the very coarse lattices with lattice spacings of up 
to 0.28 fm (corresponding to a lattice scale of 1/a ~ 0.7 GeV) this is questionable. Also, 
the range of quark masses or equivalently meson masses to which NLO-SU(2) ChPT is 
applicable will have to be determined. Eventually, the fit quality, which we measure by the 
standard x 2 /d.o.f., will be used to decide on these issues. In Fig.[3]we provide a landscape 
of the simulated meson masses at the various lattice spacings. The horizontal dashed lines 
indicate the various mass-cuts we applied in our global combined fits. One should keep in 
mind, that for the combined fit each point in the landscape plot represents two data points: 
one for the meson mass and one for the meson decay constants. 

In Fig. H] we show the result of a combined fit to the meson masses and decay constants to 
the data at all available lattice spacings and a mass range of 135 MeV < M n < 390 MeV, i.e. 
excluding only the heaviest point at each simulated /3-value. These plots show all available 
data points, whether or not included in the fit range. We marked those points which have 
been excluded by green symbols, while points included in the fit range are marked by black 
symbols. As one can already judge by eye, this fit gives a bad x 2 /d.o.f. of about 4.1, although 
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FIG. 3. Landscape plot of the simulated (squared) meson masses versus the lattice spacing a. 
The horizontal dashed lines indicate the various mass cuts we apply in our chiral fits. The blue 
dashed line indicates our final preferred fit range, see text for details. 



the fitted parameters are in a reasonable range, cf. Fig. [51 In a first step, we consider reducing 
the mass range by excluding more and more of the heavier points as indicated by the dashed 
horizontal lines in our landscape plot, Fig. [21 In Fig. Owe show how the fitted parameters and 
the x 2 /d.o.f vary when we change the mass range in the fit, still considering the ensembles 
at all available lattice spacings. First we focus on those fit ranges including the near physical 
points which are shown above the top-most dashed horizontal line in each of the plots in 
Fig-El (In the fit results shown below that line the near physical points and subsequently 
other points with light meson masses have been excluded, we will comment on those results 
below.) As one can see, the x 2 /d.o.f. improves by narrowing the fit range and for the range 
135 MeV < M n < 240 MeV it is already comparable with 1.0. Also the four fit parameters 
seem to reach a plateau at this fit range. Note, that for the two narrowest fit ranges (upper 
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FIG. 4. Combined fit for all lattice scales, 135 MeV < M w < 390 MeV. Left panel: meson decay 
constant, right panel: squared meson mass divided by quark mass ratio. Points marked by black 
symbols are included in the fit, those marked by green symbols are excluded. The physical values 
are indicated by dashed gray lines. 

mass limit at 160 or 195 MeV) the number of available data points might be too small 
as is also indicated by the larger error bars at these ranges. We already show in these 
plots the central values and error bands as determined from our preferred fit below, just to 
demonstrate that these error bands are compatible with the values obtained from fits to all 
ensembles. For completeness, in Fig.[6]we show the LECs £3, £4 and the phenomenologically 
relevant ratio f% hys /f as determined from our fits with different mass ranges. 

In a second step, we will now examine whether or not all the available lattice spacings 
already lie in the scaling region. Remember, that since the meson mass and decay constant 
define our scaling trajectory, no terms modelling discretization effects have been added to 
our chiral formulae. To test the scaling behaviour, we excluded ensembles belonging to 
one or more gauge coupling (3 from the fits. In Fig. [7] we show how the fit parameters and 
X 2 /d.o.f. change with respect to which lattice spacings are included in the fits. We show 
this for two different mass ranges, both including the near physical points. The leftmost 
point on each of these plots is from a fit to all available lattice spacings. Then, separated by 
vertical dashed lines, follow groups of fits where one, two, three or four lattice spacings have 
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FIG. 5. Results for the fitted parameters 
(top and middle panels) and x 2 /d-o.f (bot- 
tom panel) from NLO ChPT fits without 
constraints using different mass ranges but 
including all lattice spacings. The solid, 
dashed and dashed- dotted blue lines for the 
fit parameters denote the central value, sta- 
tistical and total (statistical plus systematic) 
error bands, respectively, from our preferred 
fit (in this specific set-up). 
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been excluded. The horizontal blue lines show our final estimate with error bands of the 
quantity displayed. Especially in the group where three lattice spacings have been excluded 
(the second from the right), we observe the parameters to reach a plateau by excluding the 
coarse lattice spacings. Overall, it is also reassuring that all the points fall into the error 
band of the combined statistical and systematic error, which we are going to discuss in the 
remainder now. 

From the previous discussion about excluding coarser lattice spacings from the chiral 
analysis, we decided to restrict the range in the lattice spacings to l/o > 1.6 GeV or 
a < 0.12 fm, i.e., only the ensembles at gauge couplings (3 = 3.75, 3.792, and 3.85 will 
be included. Figures |8] and M show the dependence of the fitted parameters and derived 
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FIG. 7. Fit parameters and x 2 /d.o.f. from NLO ChPT fits without constraints to two different 
mass ranges, where several lattice spacings have been excluded. Vertical dashed lines group points 
where no, one, two, three or four (from left to right) lattice spacings have been excluded. The 
solid, dashed and dashed-dotted horizontal blue lines denote the central value, statistical and total 
(statistical plus systematic) error bands, respectively, from our preferred fit. 
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quantities on the range of meson masses included in the fit range. These figures should be 
compared with Figs. |5] and |6j Again, we observe the same pattern of reaching plateaus when 
excluding more and more heavier meson masses. Eventually, we decided to take the mass 
range 135 MeV < M n < 240 MeV as our preferred fit, from which we will quote the central 
values and statistical errors. This choice is also indicated in the landscape plot, Fig. El by 
the dashed blue lines. The combined global fit has an acceptable x 2 /d.o.f. = 1.4(0.6) and is 
shown in Fig.[10j In the top panels only the data points included in the fit are plotted, while 
the bottom panels show the excluded data points as well. To estimate the systematic error 
on a fitted parameter, we take the variance of this parameter with respect to the fits using 
different mass ranges, which also include the near physical points. These are the top-most 
points above the first horizontal dashed line in each panel of Figs. El and El This procedure 
results in the set of LECs given in the left column of Tab. lIIII Note, that only the first four 
parameters are fit parameters, while the remaining ones are subsequently derived from this 
set of parameters. The central values and error bands (statistical and combined statistical 
and systematic) have always been shown in the various compilations of fit results (Figs. El 
El d El E])- It is reassuring that basically all relevant results are compatible with these error 
bands, which a posteriori justifies our procedure of estimating the systematic error. 

In the remainder of this section, we want to discuss briefly the influence of the near phys- 
ical points on the fits and the results of the fitted quantities. As an example, we provide in 
Fig.[H]the result of a combined global fit in the mass range 195 MeV < M n < 275 MeV using 
ensembles with 1 /a > 1.6 GeV. As can be seen from the insert magnifying the region around 
the physical point, for the pion decay constant the fit misses the points simulated in that re- 
gion by a couple of standard deviations resulting in an extrapolated f% hys = 129.5(0.2) MeV 
(stat. error only), which has to be compared to the value of 130.41 (0.03) (0.20) quoted by 



the PDG 



20j . For the extrapolated squared pion mass the situation looks better. In Fig. 



we provide the values for f% hjs and (M^ hys ) 2 extrapolated to the physical point by our NLO- 
ChPT fits with various mass ranges using ensembles with 1/a > 1.6 GeV only. In these 
two plots, the solid and dashed blue lines represent the central value and (total) error band 



for these quantities as quoted by [I9j, [20|. As one can clearly see, while the extrapolated 
squared pion mass is within errors compatible with the experimental result, the extrapolated 
pion decay constant starts to shift towards lower values once the nearly physical points are 
excluded from the fit range. This shift increases with both an increasing lower and higher 
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FIG. 8. Results for the fitted parameters 
(top and middle panels) and x 2 /d.o.f (bottom 
panel) from NLO ChPT fits without con- 
straints using different mass ranges and in- 
cluding lattice spacings 1/a > 1.6 GeV. The 
solid, dashed and dashed-dotted blue lines for 
the fit parameters denote the central value, 
statistical and total (statistical plus system- 
atic) error bands, respectively, from our pre- 
ferred fit. 
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FIG. 9. Results for the parameters £3, £4 (top 
panels) and the ratio f£ hys / f (bottom panel) 
from NLO ChPT fits without constraints us- 
ing different mass ranges and including lat- 
tice spacings 1/a > 1.6 GeV. The solid, 
dashed and dashed- dotted blue lines for the 
fit parameters denote the central value, sta- 
tistical and total (statistical plus systematic) 
error bands, respectively, from our preferred 
1 - 08 fit. 



bound on the fit range as can be seen from the plot. This behaviour is also in accordance 
with what could be observed in the mass range plots for the fitted parameters (Figs.0 E]) 
and derived phenomenological quantities (Figs.0 |9]) shown before. Whereas for quantities 
like ^ phys and A 3 (or £ 3 ) which predominantly influence the quark mass dependence of M%, 
the values from fits excluding ensembles with low meson masses are still in agreement with 
our estimate. This is no longer the case for the quantities / and A 4 (or £4) predominantly 
influencing the quark mass dependence of /„.. It is noteworthy though, that the ratio /*? ys // 
(where f% hys is the value extrapolated from the fit) as well shows a tendency to shift towards 
higher values once more and more lighter masses are excluded, see e.g. bottom panel of 
Fig. Mi although both / and /^ hys are shifting towards lower values, see Figs. [8] and [12l re- 
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FIG. 10. Combined fit for lattice scales 1/a > 1.6 GeV and meson masses 135 MeV < < 
240 MeV. Left panels: meson decay constant, right panels: squared meson mass divided by quark 
mass ratio. The top panels show only the points, which have been included in the fit, while the 
bottom panels show all points. There, points marked by black symbols are included in the fit, those 
marked by green symbols are excluded. The physical values are indicated by dashed gray lines. 



spectively. In our opinion, the aforementioned observations show the danger when applying 
ChPT-formulae to (lattice) data, where the lowest pion masses are still too heavy. 
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unconstrained parameter-reduced 
x ph ys/(10 -2 GeV 2 ) i. 85 78(l7) stat (39) syst 1.8639(18) stat (44) syst 



//MeV 122.70(08) stat (41) syst 122.73(06) stat (28) syst 

A 3 /MeV 628(23) stat (57) syst 678(40) stat (119) syst 

A 4 /MeV l,012(16) stat (83) syst 1, 006(15) stat (71) syst 

h 3.08(07) stat (19) syst 3.23(12) stat (30) syst 

h 4.03(03) stat (17) syst 4.02(03) stat (14) syst 

f&vsff 1.0627(05) stat (30) syst 1.0626(06) stat (24) syst 



TABLE III. Results for LECs from unconstrained (left column, see Sec. IIVB II for details) and 
parameter-reduced fits (right column, see Sec. II V B 21 for details). In case of the unconstrained 
fits the first four entries (x phyf \ f, A3, A4) are free fit parameters while the remaining entries are 
derived from these. For the parameter-reduced fit only the first two entries (x phys > /) are f ree fit 
parameters. These two sets are used to calculate the final values as quoted in Eqs. (fTM22"l) . 

2. Parameter-reduced NLO-ChPT 

As mentioned in the beginning of this section, we also considered parameter-reduced 
NLO-ChPT fits which are constrained to reproduce the physical values of the pion mass 
and decay constant at the physical quark mass, see Eqs. (FLU [T5|) . This will reduce the 
number of fit parameters in our combined global NLO-fits from four to two, leaving us with 
X phys = 2Bm phys and / as fit parameters in the constrained fits. Note that in this approach 
the low-energy constants £3 and £4 can still be determined, but they are no more independent, 
since they are linked to the LO constants as stated in equations (|T2lfT3l). The main purpose 
of this excercise is to show that our unconstrained fits did not need the additional degrees 
of freedoms to work, given that we used the pion mass and decay constants to set the 
scales for our lattice data. We performed fits with different mass ranges and limiting the 
lattice spacings included in our fits as before. In Fig. [13] we show the results for the fitted 
parameters and derived quantities as well as the x 2 /d.o.f. for different mass ranges including 
only ensembles with 1/a > 1.6 GeV. The central values and error bands shown in these 
plots have been obtained as before: the central value and statistical error is the one from 
the fit with the mass range 135 MeV < M n < 240 MeV, while the systematic error has 



23 




FIG. 11. Combined fit for lattice scales 1/a > 1.6 GeV and meson masses 195 MeV < < 
275 MeV excluding the nearly physical points. Left panels: meson decay constant, right panels: 
squared meson mass divided by quark mass ratio. The top panels show only the points, which have 
been included in the fit, while the bottom panels show all points. There, points marked by black 
symbols are included in the fit, those marked by green symbols are excluded. The physical values 
are indicated by dashed gray lines. 
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FIG. 12. Extrapolated /p^ 8 {left panel) and (MP hys ) 2 (right panel) from NLO ChPT fits without 
constraints using different mass ranges and including lattice spacings \ ja > 1.6 GeV. The solid 
and dashed blue lines denote the central value and total error bands, respectively as quoted by 

id, 



20]. 



been obtained from the variation with respect to the results from other fitting ranges in the 
meson mass. (We do not show the results from fits excluding the nearly physical points in 
this case, since due to the constraints these results now more or less agree with the results 
shown.) Figure [H] contains the plot showing the influence of the included lattice spacings on 
the fit results. We would like to state, that these parameter-reduced fits work equally well 
as the unconstrained ones, as measured by the resulting values for % 2 /d.o.f.. Our estimates 
for the LECs and derived quantities from this fitting procedure are given in the right column 
of Tab. Ill II (the first two parameters were fitted, the remaining ones were derived from these 
two). These values are in good agreement with the results from the unconstrained NLO- 
fits, given in the left column of that same table, but show a twice as large statistical and 
systematic error for A 3 or equivalently £3, whereas the remaining uncertainties are roughly 
the same or slightly reduced. 
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FIG. 13. Results for the fitted parameters (top panels) and x 2 /d.o.f (bottom right panel) from 
parameter-reduced NLO ChPT fits using different mass ranges and including lattice spacings 1/a > 
1.6 GeV. Also shown are the derived quantities £3, £4 (middle panels) and f% hys /f (bottom left 
panel). The solid, dashed and dashed-dotted blue lines for the fit parameters denote the central 
value, statistical and total (statistical plus systematic) error bands, respectively, from our preferred 
fit. 

C. Results for LECs 



As our final set of low-energy constants determined from the NLO SU(2)-ChPT fits we 
quote the following values (which supersede the preliminary results in 

X phys = 2Bm phys = 1.8609(18) stat (74) syst • 1CT 2 GeV 2 , (16) 

/ = 122.72(07) stat (35) syst MeV , (17) 

A 3 = 653(32) stat (101) syst MeV , (18) 

A 4 = 1, 009(16) stat (77) syst MeV , (19) 
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FIG. 14. Fit parameters and x 2 /d.o.f. from parameter-reduced NLO ChPT fits to two different 
mass ranges, where several lattice spacings have been excluded. Vertical dashed lines group points 
where no, one, two, three or four (from left to right) lattice spacings have been excluded. The 
solid, dashed and dashed-dotted horizontal blue lines denote the central value, statistical and total 
(statistical plus systematic) error bands, respectively, from our preferred fit. 
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% = 3.16(10) stat (29) syst j 
? 4 = 4.03(03) 8t at(16) syst j 



/f y 7/ = 1.0627(06) stat (27) syst . 



(20) 
(21) 
(22) 



Here we averaged the central values and statistical uncertainties from the unconstrained 
(Sec. II VB lj) and parameter-reduced (Sec. IIVB 2j) fits (both summarized in Tab. IHIj) . For 
the square of the systematic uncertainty we sum the squares of the average systematic 
uncertainty and of the spread of the central values. 



V. FITS TO NNLO-SU(2) CHPT 

As a further step, we examined the effects of fitting our data to SU(2)-ChPT including 
terms up to next-to-next-to- leading order (NNLO). The motivation for this excercise is two- 
fold. On the one hand, we would like to check whether or not our results obtained for the 
LECs by fitting to NLO-ChPT are stable when considering the next higher order in ChPT. 
On the other hand, there is some general interest in whether the amount of data is sufficient 
to reliably determine the additional fit parameters appearing at NNLO and whether NNLO- 
ChPT is superior to NLO-ChPT in terms of describing the data (as measured, e.g., by the 
X 2 /d.o.f.). 

A. Fit formulae 



The formulae for the squared meson mass and the meson decay constant up to NNLO 
read (cf. [la, ba|) 



u = s 

NNLO m 2 = ( 



NNLO 



1 + 
1 - 

Xi 
16n 2 f 2 

Xi 
16tt 2 / 2 



X log^r + NNLO M 2 



16vr 2 / 2 ~° A 2 



X 



X 



^ P l0g M + NNLO, 



1 

306 



— j 60 log ^- - 91og^ - 49 J + Ak M 2 



X 



(23) 
(24) 

(25) 



X 



X 



X 



180 



301og^- + 61og-^ - 61og^ - 231 +4% 



(26) 
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Up to this order three new fit parameters enter: one additional low energy scale A 12 and 
the two NNLO-LECs k M 2, kf. The low energy scale A12 is related to the low energy scales 
usually denoted by Ai and A 2 in the literature via 

logA? 2 = ^logA? + A i g a|. (27) 

The low energy scales Ai, A 2 already appear separately in the NLO-formulae for other 
quantities (e.g. scattering lengths in 7T7r-scattering but since in our case only the 

combination Ai 2 appears, we will not be able to distinguish between them. As before, the 
low-energy scales can also be expressed via the LECs li as in Eq. ([TP]) , i.e. 

*< = lo S^^' , = 1,2,3,4,12. (28) 



B. Combined global fits 

We applied the same fit strategy for the NNLO-fits as for the NLO-fits (see previous 
section): the scale and light quark mass at the physical point as determined in Sec. lIIII 
are used and combined global fits to our data for the meson masses and decay constants 
are performed. First, we will discuss fits without any assumptions on the fit parameters 
(Sec. lVBTT) . Later we will also constrain the additional fit parameters entering at NNLO by 
phenomenologically motivated estimates (Sec. IV B"2l) . 



1. NNLO fits without priors 

We will start our discussion with fit ranges including the nearly physical points and only 
consider varying the upper mass limit of the fit range. Since now three more parameters 
have to be determined, it can be expected that we will have to include more data-points, 
i.e., including higher meson masses, compared to the NLO-case. Indeed, with too few data- 
points either the fitter could not find a solution at all or some of the fitted parameters had 
big numerical uncertainties. This can be seen in the compilation in Fig.HS] of fit results 
using ensembles with 1/a > 1.6 GeV, e.g., for the fit ranges 135 MeV < M n < 240 MeV, 
135 MeV < M n < 275 MeV and maybe also 135 MeV < M n < 340 MeV. 

The plots of a sample fit with the range 135 MeV < M n < 390 MeV, 1/a > 1.6 GeV 
are given in Fig. [16] (only data points included in the fit range are shown in these plots). 
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FIG. 15. Results for the fitted parameters and x 2 /d.o.f from NNLO ChPT fits without constraints 
using different mass ranges and including lattice spacings 1/a > 1.6 GeV. The solid, dashed 
and dashed- dotted blue lines for the fit parameters denote the central value, statistical and total 
(statistical plus systematic) error bands, respectively, from our NLO-fits. 
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FIG. 16. Combined NNLO-fit for lattice scales 1/a > 1.6 GeV and meson masses 135 MeV < 
M n < 390 MeV. Left panel: meson decay constant, right panel: squared meson mass divided by 
quark mass ratio. The solid blue lines show the complete (up to NNLO) fit, whereas the dashed 
red lines show the LO+NLO contribution of the full NNLO-fit. Only data points included in the 
fit range are depicted in the plots, the physical values are marked by dashed gray lines. 

Although technically the fit seems to work well resulting in an acceptable x 2 /d.o.f. = 1.8(0.6) 
one should take the result cum grano salis. Besides the fitted curve up to NNLO (solid blue 
line with error band indicated by blue dotted lines), each plot also shows the LO+NLO 
contribution separately (dashed red line with error band indicated by red dotted lines), 
i.e., the NNLO contribution being the difference between the solid blue and the dashed red 
line. In this case, the NLO- and NNLO-contributions taken separately seem to have big 
uncertainties. In other cases, we found the NNLO-contribution to have an unnaturally large 
effect even at small quark masses. 

At this point, we have to conclude that NNLO-fits to our current data are not convincing 
for the reasons outlined above. This situation might improve, once more data points in the 
region between the physical point (or below) and, say, 200 MeV or 250 MeV will be added 
to the analysis, allowing for fits using a smaller mass range (like we were able to do for the 
NLO-fits). Since the generation of such data points is not planned for the near future, we 
will in the remainder of this section examine whether constraining the additional NNLO fit 
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parameters can serve as a remedy to this situation (and which side-effects this remedy has). 

2. NNLO fits using priors 

To stabilize the NNLO fits, we examined the effect of using priors for the three additional 
fit parameters A 12 , &m 2 > and kf. Since we are mainly interested to learn how the parameters 
already appearing at NLO change, when going from NLO to NNLO fits, we did not consider 
adding priors for these parameters as well. 

The priors were implemented in the fitting procedure via an augmented x 2 - Instead 
of minimizing the usual \ 2 (sum of deviations of the data points from the fitted function 
weighted by the uncertainty of that data point) the modified 

2 2 (Pn ~ p pvioT ) 2 

^augmented X ^ , / prior\2 (^^) 

n \Gn ) 

is minimized by the fitting routine. Here n runs over all fit parameters p n to which a 
prior pP rior with width u^ lov has been assigned. That is, the fitting routine "punishes" an 
parameter for leaving the prior intervall p^ rior ± cjP rior . 

In the following we will use either a prior for £12 alone or together with priors for ^2 
and kf. We will assign the following values for the priors to the fit parameters: 

£?f r = 2.1 ± 0.3 (30) 
k^T = ± 10 (31) 
kf OT = ± 10 (32) 

The prior for £ 12 has been obtained from the estimates for the LECs i\ = —0.4 ± 0.6 and 
I2 = 4.3 ± 0.1 extracted from 7T7r-scattering data in {25] via 

7 H 

/jprior ' «prior . /jprior 

cf. Eqs. (1271) and (!28l) . This would translate into a prior for the low energy scale 

Ag ior = 385 MeV ± 58MeV. (33) 

The priors for the NNLO-LECs k M i and kf are merely based on assuming a natural order 
of magnitude for these corrections. We are aware that the latter is a rather weak argument, 
for that reason we did not use the priors on fc M 2 and kf alone and results from fits using 
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FIG. 17. Results for the fitted parameters and x 2 /d.o.f from NNLO ChPT fits with a prior for 
A12 using different mass ranges and including lattice spacings 1/a > 1.6 GeV. The solid, dashed 
and dashed- dotted blue lines for the fit parameters denote the central value, statistical and total 
(statistical plus systematic) error bands, respectively, from our NLO-fits. The used prior and its 
width is indicated by the dashed black line and the shaded grey area, respectively. 
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FIG. 18. Results for the fitted parameters and x 2 /d.o.f from NNLO ChPT fits with priors for A12, 
kj^2, and kf using different mass ranges and including lattice spacings l/o > 1.6 GeV. The solid, 
dashed and dashed-dotted blue lines for the fit parameters denote the central value, statistical and 
total (statistical plus systematic) error bands, respectively, from our NLO-fits. The used priors 
and their widths are indicated by the dashed black lines and the shaded grey areas, respectively. 
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FIG. 19. Combined NNLO-fit for lattice scales 1/a > 1.6 GeV and meson masses 135 MeV < 
M w < 340 MeV using a prior for A12. Left panel: meson decay constant, right panel: squared 
meson mass divided by quark mass ratio. The solid blue lines show the complete (up to NNLO) 
fit, whereas the dashed red lines show the LO+NLO contribution of the full NNLO-fit. Only data 
points included in the fit range are depicted in the plots, the physical values are marked by dashed 
gray lines. 

these priors should be taken with some caution. However, these priors are to some extent 
justified by the NNLO fits where no priors have been used (see Fig.[15]) and also, as we will 
see, by fits where only the prior on £ 12 has been used. 

The results from the fits to different mass ranges using priors on £ 12 only or on /c M 2 
and kf as well are shown in Figs. [T7] and HB] respectively, where only the ensembles with 
1/a > 1.6 GeV have been used. The used priors are indicated by the gray shaded areas. 
As an effect, now fits to smaller mass ranges are possible and/or are more stable judging 
from the uncertainties of the fit parameters. In Figs. [19] and [20] we show as examples the 
fits for the mass range 135 MeV < < 340 MeV. Whereas in both cases for the decay 
constant the NNLO-contribution seems to be reasonable small, for the meson mass the 
NNLO-contribution is more substantial. Nevertheless, the error bands for the LO+NLO 
contribution are now reasonable as well, showing again the stabilizing effect of using the 
priors. 
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FIG. 20. Combined NNLO-fit for lattice scales 1/a > 1.6 GeV and meson masses 135 MeV < 
M w < 340 MeV using priors for A12, k M 2, and kf. Left panel: meson decay constant, right panel: 
squared meson mass divided by quark mass ratio. The solid blue lines show the complete (up to 
NNLO) fit, whereas the dashed red lines show the LO+NLO contribution of the full NNLO-fit. 
Only data points included in the fit range are depicted in the plots, the physical values are marked 
by dashed gray lines. 



In Figs. 1211 and 1221 we compare the LECs and the ratio f% hys /f from the NNLO fits to 
different mass ranges using either no priors or one of our choices for the set of priors (all for 
ensembles with 1/a > 1.6 GeV). Our final results with total error bands from the NLO-fits 
are always indicated by the blue lines. In the case of £3, the fits using priors shift the value 
of the LECs a bit upward, whereas £4 and f% hys /f just fluctuate within the total error band 
from the NLO-fit. The comparison for the LEC £12 again demonstrates the stabilizing effect 
of the priors or — in the other words — the difficulties we encountered in the fits without using 
priors. In the latter case (top right panel of Fig. l22l) . the fit result only comes close to the 
phenomenological estimate, Eq. ( 130|) . for fit ranges including meson masses of 390 MeV or 
higher. 
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FIG. 21. Results for the LECs £3 (left panels) and ^4 (right panels) from NNLO fits to ensembles 
with 1/a > 1.6 GeV and different mass ranges. The top, middle, and bottom panels show the 
results using no priors, a prior for A12, and priors for A12, and kf, respectively. The solid, 

dashed, and dashed dotted blue lines indicate the central value, statistical and total (statistical plus 
systematic) error bands from our NLO-fits. 



VI. CONCLUSIONS 



In this paper we presented a determination of the NLO low-energy constants £3 and £4 of 
SU(2) chiral perturbation theory from 2+1 flavor lattice simulations with staggered fermions. 
In addition we gave results for the LO quantities / (or f% hys /f) and B (or x phys = Bm phys ) . 
The quantities £ 3j4 are also expressed in terms of the scales A 3 4 , respectively. 
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FIG. 22. Results for the ratio /P hys // (left panels) and the LEC £ 12 (right panels) from NNLO fits to 
ensembles with 1/a > 1.6 GeV and different mass ranges. The top, middle, and bottom panels show 
the results using no priors, a prior for A12, and priors for A12, &m 2 > an d kf, respectively. The solid, 
dashed, and dashed dotted blue lines indicate the central value, statistical and total (statistical plus 
systematic) error bands from our NLO-fits. The dashed black lines and grey shaded areas indicate 
the used prior and its width (where applicable). 



Our final results as presented in equations f|T6l - [22l) stem from fits which use the NLO 
functional form. The result for ^ ph y s amounts to a condensate parameter 

B(2 GeV) = 2.682(36) (39) GeV (34) 



if one divides out the value of the average light quark mass from Ref . 26j, |27j . Moreover, 
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after multiplying with F 2 = f 2 /2 from (fT7|) one obtains 

E(2 GeV) = 2.020(27)(31) GeV 3 or E(2 GeV) 1/3 = 0.2723(12) (14) GeV (35) 
where all quantities given at the scale /i = 2 GeV refer to the MS scheme. These results are 



reasonably consistent with the high-quality entries in Table 10 of Ref. [19], in particular with 
those by the MILC [17], RBC/UKQCD (28f] , and ETM I29J collaborations, to mention some 
of the most precise determinations. Similarly, our results 

4 = 3.16(10) (29) and l A = 4.03(03)(16) 

as stated in (120]) and (l2Tj) . respectively, are in good agreement with the broad majority of 



;he entries in Table 12 of Ref. 
and ETM 



19|, in particular with those by the MILC |30|, RBC/UKQCD 



29| collaborations. Finally, our result fT2"2"|) for the ratio f% hys / f agrees with 



the entries of Table 1 1 in of Ref. 



19] , perhaps with some slight tension when compared to 



3l|. 



the recent 2+1+1 flavor result by the ETM collaboration 

We have carefully examined the effect of various cuts on the data, in particular the effect 
of requesting a -1 > 1.6 GeV and the effect of limiting the pion mass range that enters the 
chiral fit. It turns out that the restriction to fine lattices improves the quality of the fits, 
and with this restriction reasonably stable NLO results are obtained for pion mass windows 
up to ~400 MeV. Given the fine grained set of pion masses available in our ensemble basis, 
we can even explore the effect of dropping some of the lighter data points. We find that £3 
is much more robust in this respect than £4 (or f% hys /f), as is evident from Fig.|9j 

Finally, we have explored the effect of adding the NNLO contribution to the functional 
ansatz. To prevent a dramatic increase of the number of free parameters we add priors 
for the new NLO combinations (in which we are not interested) and the genuine NNLO 
coefficients (to which our data show little sensitivity). We find it very reassuring that these 
prior-aided NNLO fits remain stable (for a reasonable range of lattice spacings and pion mass 
windows), and that the resulting fits show a very natural ordering between LO, NLO, and 
NNLO contributions (out to ~ 400 MeV) . Moreover, the NLO coefficients I 3 and I A as 
determined from these prior-aided NNLO fits are in good agreement with the results of the 
direct NLO fits. We take this as a sign that our assessment of the systematic uncertainties 
of these quantities is true and fair. 

As stated in SectionQV] we did not attempt to include terms designed to absorb cut-off 
effects into our formulas. The purpose of the present work was to explore whether such 
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terms are mandatory; we find that they are not (with the level of precision of our data), 
albeit at the price of pruning the dataset to include only lattices with a -1 > 1.6 GeV. It 
remains open to which extent it would prove helpful to include such cut-off terms; we may 
address this point in a future study. 

In summary, our results (fl6l- l22l) and (134 11351) indicate that SU(2) chiral low-energy con- 
stants can be determined on the lattice with a precision at the level of a few percent for the 
LO quantities B, E, /, and the level of 0(10%) for the NLO scales A 3 , A 4 . 
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